Applying a new method to the 2-d Ising transition 
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The Ising ferromagnetic model on a square lattice is revisited using the Galam Unifying Frame 
(GUF), set to investigate the dynamics of two-state variable systems within the frame of opinion 
dynamics. When combined with Metropolis dynamics, an unexpected intermediate "dis/order" 
phase is found with the coexistence of two attractors associated respectively to an ordered and a 
disordered phases. The basin of attraction of initial conditions for the disordered phase attractor 
starts from zero size at a first critical temperature T c i to embody the total landscape of initial 
conditions at a second critical temperature T C 2 with T c i « 1.59 and T C 2 « 2.11 in J/fcs units. It 
appears that T C 2 is close to the Onsager result T c ~ 2.27. The transition, which is first-order like, 
exhibits a vertical jump to the disorder phase at T C 2, reminiscent of the rather abrupt vanishing of the 
corresponding Onsager second order transition. However, using Glauber dynamics combined with 
GUF does not yield the intermediate phase and instead the expected single transition at T c « 3.09. 
Accordingly, although the "dis/order" phase produced by the GUF - Metropolis combination is 
not physical, it is an intriguing result to be understood. In particular the fact that Glauber and 
Metropolis dynamics yield so different results using GUF needs an explanation. The possibility of 
extending GUF to larger clusters is discussed. 



I. INTRODUCTION 

The Ising model is a seminal model of statistical 
physics, which has inspired thousands of scientific papers 
[!l,|2| . In particular the two-dimensional nearest neighbor 
square lattice ferromagnetic Ising model is a cornerstone 
to the understanding of collective phenomena being the 
unique case to exhibit a second order phase transition 
with an exact analytical solution ||. All its properties 
are considered to be known. 

Accordingly, in this work we use the two-dimensional 
nearest neighbor ferromagnetic Ising model to investi- 
gate the properties of the Galam Unifying Frame (GUF), 
set to investigate the dynamics of two-state variable sys- 
tems within the frame of opinion dynamics [|[ . Applying 
GUF combined with Metropolis dynamics for the prob- 
ability of a single spin flip [5|, Q , reveals an unexpected 
twofold order-disorder Ising transition. The transition is 
no longer at once from the ordered state into the dis- 
ordered state but instead goes through an intermediate 
phase denoted "dis/order" phase, for which two attrac- 
tors exist, each one being associated to a different basin 
of attraction of initial conditions. Their respective sizes 
are a function of the temperature. 

The basin of initial conditions to the disordered phase 
starts from zero size at a first critical temperature T c \ sa 
1.59 to end up embodying the total landscape of initial 
conditions at a second critical temperature T C 2 ~ 2.11. 
Temperature values are given in J/fcs units where J is 
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the positive coupling constant and ks the Boltzmami 
constant. To make the presentation lighter from now 
on we take J/ks = 1- 

For the Ising system, noting p the proportion of spins 
in state +1 and (1 — p) the proportion of spins in state 
— 1, three different regimes are obtained for the dynamics 
as a function of temperature. 

(i) Starting at T = the phase diagram contains two 
attractors p_ = and p + = 1 with a separator p c = \. 
Increasing temperature moves the two attractors toward 
the separator with p_ > and p+ < 1 . 

(ii) Then, at a temperature T c \ the separator p c = \ 
turns to an attractor. All these steps are expected from 
the usual description of the order/disorder Ising transi- 
tion. However, while in the classical scheme, the change 
of status from separator to attractor is tuned on by the 
simultaneous merging of the two attractors p_ and p + at 
p c = ^, here the opposite process occurs. Two new sym- 
metrical separators p c _ < | and p c+ = 1 — p c _ > i are 
expelled from p c = which becomes a third attractor. 
That is an unexpected process. 

(iii) Keeping on increasing the temperature, on each 
side of the attractor p c = |, the corresponding attractor 
and the new separator move towards each other to even- 
tually coalesce at a second critical temperature T c i and 
then disappear. In the range T > T C 2 only the attractor 
Pc = I exists. 

At this stage it happens that T c2 ~ 2.11 is close to 
the Onsager exact result T c ss 2.27. Moreover, the tran- 
sition, which is first-order like, exhibits a vertical jump 
to the disorder phase at T C 2, reminiscent of the rather 
abrupt vanishing of the corresponding Onsager second 
order transition. 

However performing the same calculations as above us- 
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ing combining GUF with Glauber dynamics instead of 
Metropolis, restores the expected phase diagram of one 
single transition from the ordered phase into the disor- 
dered phase. The transition is continuous like and occurs 
at T c rj 3.09 against T c = 4 for the Mean Field counter- 
part. 

Accordingly, although the "dis/order" phase produced 
by the GUF - Metropolis combination is not physical, it 
is an intriguing result to be understood. In particular 
the fact that Glauber and Metropolis dynamics yield so 
different results using GUF needs an explanation. The 
possibility of extending GUF to larger clusters is also 
briefly discussed. 



II. APPLYING THE GUF DYNAMICS 

With no external field, the Hamiltonian for the square 
lattice ferromagnetic Ising system is given by 



H 



M 



(1) 



a sum over all nearest neighbors i and j 
= ±1. The associated magnetization is M = 
<7i = 2p — 1 where N is the total number of spins. 



where [ij] 
and <7j = 

N i 

The problem has been solved exactly with T c = 
2/arcsinh(l) w 2.27. In contrast, the classical mean field 
approach yields T c = 4. However in higher dimensions 
or with an external field, Monte Carlo (MC) simulations 
are required. They are implemented using a dynamics 
which respects detailed balance. The most common ones 
arc the Glauber and Metropolis dynamics [!, [f| . 

In the Glauber dynamics, all spins are investigated in 
a sequential manner for each MC step. Given a spin cr^ 
in a configuration /Zj with energy , the flip cr.; — > — er, 
creates the new configuration r)i with energy E, h = —E fii . 
The actual flip is implemented with the probability 



1 



l + cxp(-2^ i /fc B T) 



(2) 



In Metropolis scheme, the spin is selected randomly 
and a MC step corresponds to N updates. The flip prob- 
ability to the new configuration is given by 



the update. Here 5 spins are used, the central one plus 
its 4 nearest neighbors yielding 



p(*+i) = X>*[p(*)]*[i-p(*)] 



5— fe 



(4) 



k=0 



where p{t) is the proportion of +1 spins at time t and 
p(t + 1) is the proportion of +1 spins at time t + 1 
after the equivalent of one MC step. The product 
p(t)] k [l — p(t)] 5 ~ k is the probability to have a group 
of five spins with k spins +1 and (5 — k) spins — 1, 
while the coefficient gk is the probability that the con- 
figuration ends up with the middle spin in the state +1 
with the other 4 unchanged following either Glauber of 
Metropolis rule. By up-down symmetry g$ = \ — g$, 
gi = 5 — 54 and (72 = 10 — 33. Performing the enumer- 
ation of the various configurations leads to respectively 
g 5 = 1 - b 2 ; g 4 = 5 - 46; g 3 = 4 and g 5 = 1/(1 + 6 2 ); g A = 
1/(1 + b 2 ) + 4/(1 + b);g 3 = 3 + 1/(1 + b) for respec- 
tive Metropolis and Glauber rules where temperature is 
incorporated using the parameter b = e~ 4 / T . Plugging 
those coefficients into Eq. (UJ) leads to respectively 

Pm = P 5 (l - b 2 ) +p\l - p)(5 - 46) + 4p 3 (l - p) 2 
+ 6p 2 (l - pf + 46p(l - P ) 4 + 6 2 (1 - p) 5 . (5) 



and 

PG = P 



3n , 4 )+p 2 (1 _ p) 3 {7 _ 4 



(G) 



+ p(l- P ) 4 (5 



1 + 6 7 
1 



1 + 6 2 1 + 6 



) + (l-p) 5 (l 



1 + 6 2 " 



M^, H = min{l,cxp[2^/fc B T]}. 



(3) 



where p — pit) with pm = pit + 1) using the Metropoilis 
rule and pa = p(t + 1) using the Glauber rule. 



III. UNCOVERING A TWOFOLD TRANSITION 

To investigate the dynamics and phase transitions pro- 
duced by Eqs. ([5j [6]) we need to solve the fixed point 
equation p{t + 1) = p(t) for each case. The attractors 
correspond to the equilibrium states while the separa- 
tor determines the flow direction when implementing the 
dynamical update. 



In parallel a new scheme has been developed in the re- 
cent years to unify the large spectrum of discrete models 
proposed to study opinion dynamics. Indeed the Galam 
Unifying Scheme (GUF) was shown to embody all avail- 
able discrete models of opinion dynamics Q. We now 
apply it to the classical ferromagnetic Ising model on a 
square lattice. 

The GUF considers a polynomial development whose 
degree is given by the size of the group used to perform 



A. GUF/Metropilis 

From pm = P we get always the solution p c = 1/2 with 
3 different regimes. One with p c = 1/2 being the unique 
attractor. The second and third ones have respectively 
two and four symmetrical solutions up and bellow p = 
1/2 as shown in Figure [TJ In the second case p c = 1/2 
is a separator while it is an attractor in the third case. 
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separator. However, they do not reach the zero magne- 
tization as b — > b C2 . Instead, they meet another pair of 
solutions 



< T < r, 



T = T, 
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FIG. 1. Landscape of the attractors and separators as a func- 
tion of temperature. The first top corresponds to T = with 
two attractors p=0, p+ = 1 and a separator p c — 1/2. Sec- 
ond from top corresponds to < T < T c \ ~ 1.59: same as 
before with p_ > and p+ < 1. Third from top: at T = T c i 
the separator is about to turn into an attractor giving rise to 
two symmetrical moving separators on each side. Third from 
bottom: for T c \ < T < T C 2 three attractors p_, p+, p c sepa- 
rated by two the separators p c -, Pc+- Second from bottom: 
at T = T C 2 the two separators p c _ and p c + coalesce with the 
two attractors p_, p+ to suppress them. Bottom: only one 
attractor at p c = 1/2 for T > T C 2 ~ 2.11, which is close to 
the exact result 2.27. 



More precisely, above 



bc 2 = \{ 1 - 16 ( 
0.150, 



3(5^/249- 9) y 3 7 v 3 7 v ' 



(7) 



we observe only the zero-magnctization solution p = 1/2. 
However, below b C2 , in the region defined by < b < 



b Cl with 



(v 41 — 6) ps 0.081, 



we have two solutions for p given by 



P± = h { C ±VC' 2 -4C(D-E) 



(8) 



(9) 



where C = (2b 2 - 8b + 6), D = {3b 2 -46+1) and 
E = V5& 4 -86 3 + 106 2 - 86+l. These solutions are at- 
tractors. They start from the values p_ = and p+ = 1 
at b = and move towards p c = 1/2 which behaves as a 



= 2c\C±y/C 2 -AC{D + E), 



that exist only in the interval 



b Cl <b<b C2 . 



(10) 



(11) 



As soon those two solutions p c ± appear, they behave 
as separators and move towards p c = 1/2, which turns to 



an attractor as b 



The values b Cl and b C2 yield the 



critical temperatures of T Cl ps 1.59 and T C2 ss 2.11 using 
T = -4/ In b. 

We see that as soon as b < b C2 , p c = 1 /2 is no longer the 
only solution, leading us to believe that this is the actual 
phase transition. This provides a much better estimate 
of the critical temperature than that of mean field theory, 
since the correct value is T c ps 2.27. Notice that T C1 ps 
1.59 had been previously observed for another attempt 
of estimating the critical temperature in the context of 
the GUF framework 0] , however missing the existence of 
T C2 . 

The series of diagrams of Figure [1] shows the moving 
landscape of attractors and separators as a function of 
varying the temperature. Two attractors p_ = and 
p + = 1 with a separator p c = | are found at T = 
as expected. Increasing temperature the two attractors 
move toward the separator with p_ > and p + < 1 also 
as expected. However at a temperature T c \ ps 1.59 the 
separator turns to an attractor p c = \ with the simul- 
taneous appearance of two new symmetrical separators 
p c _ < i and p c+ = 1 — p c - > \- That is an unexpected 
result. 

Increasing still the temperature, on each side of the 
attractor p c = h, the attractor and the new separator 
move towards each other to eventually coalesce at T cl ps 
2.11 and then disappear. At T > T C 2 only the attractor 
p c = 2 exists. 

Figure [2] shows the variation of all the five fixed points 
as a function of b = cxp~~. The arrows shows the flow 
direction while iterating the local updates. Dark solid 
lines correspond to attractors, i.e., p+,p~ and p c for b > 
b c i- Dashed lines are separators, i.e., p c +,p c - and p c for 
b < b c i 

Figure [3] shows the three kind of dynamics obtained to 
describe the GUF twofold Ising transition. At b = 0.006 
and b = 0.16 stand the classical typologies with all initial 
configurations ending up at at an ordered phase with 
either a positive (po > 0.50) or negative (po < 0.50) 
magnetization for the first case and to zero magnetization 
for the second case. The Figures at b = 0.10 and b = 
0.13 are unexpected with two possible equilibrium states 
depending on the initial configuration. It is worth to 
stress that the equilibrium state is not probabilistic. It 
is either M ^ or M = depending on the value of Pq. 
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FIG. 2. The variation of all the five GUF fixed points as a 
function of b = exp(—4:/T) using Metropolis dynamics. The 
arrows show the flow direction while iterating the local up- 
dates. Dark solid lines correspond to attractors, i.e., p+,p- 
and pc for b > b c i- Dashed lines are separators, i.e., p c +,p c _ 
and p c for b < b c i. 






FIG. 3. The y axis shows the variation of po under repeated 
application of the update rule pg using Eq. ([5]) from po = Of 
to po — 0.90 with an incrustation of +0.025. The x axis shows 
the number of iterations up to 80. The series are shown for 
b = 0.006, O.fO, 0.13, 0.16. Up left (6 = 0.006) corresponds to 
the usual ordered phase while bottom right (6 = 0.f6) to the 
usual disordered phase. Upright (6 = O.fO) and bottom left 
(b — 0.f3) show the uncovered new transition. The number 
of iterations are labeled along the abscise with 0, 20, 40, 60, 



B. GUF/Glauber 

We now repeat above calculations using the Glauber 
scheme (Eq. [2]) instead of Metropolis (Eq. [3]). From Eq. 
[5] solving pc = p yields only three solutions, p c = 1/2 
and 



b)iVT 



(12) 



were S = -1 + 36 - 3b 2 + b 3 and T = -1 + 3b + b 2 + 5b 3 . 
The separators p c ± obtained in the precedent Subsection 



FIG. 4. Landscape of the attractors and separators as a func- 
tion of temperature. The first top corresponds to T = with 
two attractors p=0, p+ = 1 and a separator p c = 1/2. Second 
from top corresponds to < T < T c ~ 3.09: same as before 
with p- > and p+ < 1. Third from top: at T = T c , p- and 
p+ merge with the separator p c , which in turn becomes an 
attractor. Bottom: only one attractor at p c = 1/2 for T > T c 



do not exist any longer making the phase diagram shown 
in Figure [4] very different from the one from Figure [T] 
However it is more like what is expected for the two- 
dimensional nearest neighbor ferromagnetic Ising model 
with no "dis/ordcr" phase and a single transition from 
the ordered phase into the disordered. 

The associated critical temperature is obtained when 
p + = P- = p c = 1/2 yielding b c ss 0.274, which gives 
T c rj 3.09, whose value is in between the mean field result 
T r — 4 and the exact result T r w 2.27. 



C. GUF/Glauber/Onsager/MF 

To compare the results obtained with respectively 
Metropolis and Glauber, the associated attractor and 
separator curves are exhibited in Figure [5] as a function 
of temperature via the variable b in abscise together with 
the Onsager exact solution for the magnetization Q 

mo = (l-sinh(2/Tr 4 )^, (13) 
where 2 /T = -1/(2 In b) and the Mean Field formula 

( 4wa/f ' 



tumf = tanh 



V T 



(14) 



using p = (m + l)/2 with M denoting Metropolis, G 
Glauber, O Onsager and MF Mean Field. Five critical 
values b c are obtained with (1, 2) for Metropolis, (3) for 
Onsager, (4) for Glauber and (5) for Mean Field. 



FIG. 5. Landscape of the attractors and separators for p as 
a function of temperature via the variable b = exp(— 4/T) in 
the abscise with Metropolis (M), Glauber (G), Onsager (O) 
and Mean Field (MF). Solid lines represents attractors while 
dashed lines are separators. However the line p c — 1/2 is a 
separator for M, G, O and MF only when b < b c i (1). For 
fcci(l) < b < 6 2(2), p c is an attractor for M and a separator 
for O, G and MF. However, for M other attractors coexist. 
When 6 C 2(2) < b < 6 C (3), p c is the unique attractor for M 
and still a separator for O, G and MF. For 6 C (3) < 6 < 6 C (4), 
p c is an attractor for M, 0, but still a separator for G, MF. 
When 6 C (4) < b < 6 C (5), p c is an attractor for M, 0, G, 
but still a separator for MF. Only when b > 6 C (5), p c = 
1/2 is the unique attractor for M, O, G and MF. We have 
bd ~ 0.081, 6 c2 w 0.150, 6 C (3) ps 0.171, 6 C (4) w 0.274, 6 C (5) ps 
0.368 yielding respectively T c i ~ 1.59, T c2 w 2.11, T c (3) « 
2.27, T c (4) « 3.09, T c (5) = 4. 

From Figure [5] it is worth to notice that Metropolis 
yields a first-order like transition with a critical temper- 
ature b C 2 ps 0.150, which is close to Onsager exact re- 
sult 6 C (3) ps 0.171, in addition to exhibit a vertical jump 
to the disorder phase at 6 C 2, reminiscent of the rather 
abrupt vanishing of the corresponding Onsager second 
order transition. On the other hand, Glauber is more 
like MF with an improvement in the value of the critical 
temperature at 6 C (4) ~ 0.274 instead of 6 C (5) ps 0.368 for 
Mean Field. 

Solid lines represents attractors while dashed lines are 
separators. However, the line p c = 1/2 is a separa- 
tor for M, G, O and MF only when 6 < 6 cl (1). For 
6 c i(l) < 6 < 6 c2 (2), p c is an attractor for M and a sep- 
arator for O, G and MF. But for M other attractors co- 
exist. When 6 c2 (2) < b < 6 C (3), p c is the unique attrac- 
tor for M and still a separator for O, G and MF. For 
6 C (3) < b < 6 C (4), p c is an attractor for M, O, but still 
a separator for G, MF. When 6 C (4) < b < b c (5), p c is 
an attractor for M, O, G, but still a separator for MF. 
Only when b > b c (5), p c = 1/2 is the unique attrac- 
tor for M, O, G and MF. We have b c i ps 0.081, 6 c2 ps 
0.150, 6 C (3) ss 0.171, 6 C (4) ps 0.274, 6 C (5) ps 0.368 yielding 
respectively T cl ps 1.59, T c2 ps 2.11, T c (3) ps 2.27, T c (4) ps 
3.09, T c (5) = 4. 
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FIG. 6. The variation of all the five GUF fixed points as 

_ _4 

a function of b — exp t . The arrows show the flow direc- 
tion while iterating the local updates. G denotes the wrong 
Glauber scheme while M is for Metroplolis. 

D. GUF/Glauber bis 

In an earlier version of the paper, a misprint in Eq. 
© in the coefficient g 3 = (3 + j^) of p 3 (l - p) 2 had 
been carried on during the iteration of the calculations 
previously done with Metropolis, leading to a very similar 
phase diagram as obtained with Metropolis in Figure [6l 

More precisely, putting jj^ instead of in (73 , which 
implies the same change in g 2 = 10 — (73 recovers the 
twofold phase diagram one obtained using Metropolis 
with 



(-7 61 U ' — - + (11(1147+ Vl059)) 1/3 ^) 

V 1147 + J 

ps 0.079, (15) 

and b C2 is the solution of 

- 1 - 66 + 3b 2 - 32b 3 + 3& 4 - 186 5 + b 6 = 0, (16) 

which cannot be solved analytically. However, the nu- 
merical solution is b C2 = 0.158. With respect to the four 
fixed points, we have 

P± = ^( F± VF 2 -mD-H)) (17) 

where F = (-6b 2 - 6), G = (b 3 - 3b 2 + 36 - 1) and 
H = Vb 6 - 18b 5 + 36 4 - 326 3 + 3b 2 -66+1. 

As for Metropolis these solutions are attractors start- 
ing from the values p_ = and p + = 1 at 6 = and 
moving towards p c = 1/2 which behaves as a separator 
as 6 is increased. At 6 = 6 C2 they coalesce with 

Pc± = ^p(F± ^F 2 -AF{G + H)) . (18) 
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The critical temperatures associated to b Cl and b C2 are 
T C1 as 1.57 and T C2 as 2.17, which is rather close to the 
exact value T C2 as 2.27. 

Although those results are wrong since resulting from 
a missing factor in Eq. ^ , they are worth to be noticed 
due to the surprising identity with the ones obtained us- 
ing Metropolis as illustrated in Figure [5] Such a coinci- 
dence is worth more investigation. Indeed, the finding of 
the "dis/order" phase was unexpected using Metropolis 
thus prompting the calculations to be carefully checked 
again and again from every part. However, the fact that 
Glauber was found to yield the same unexpected result 
as Metropolis, was expected, and thus made the misprint 
unnoticed in the coefficient 173 in Eq. (jB]) with a recheck- 
ing of only the calculations starting from Eq. (JB]). 



IV. DISCUSSION 

In this work, the application of the GUF combined 
with either Metropolis or Glauber schemes to the clas- 
sical two-dimensional ferromagnetic Ising model, was 
shown to exhibit some rather unexpected results. In par- 
ticular the Metropolis scheme leads to the appearance of 
an intermediate "dis/order" phase between the ordered 
and disordered phases, turning first-order like the associ- 
ated transition. It happens that the corresponding crit- 
ical temperature T c i ~ 2.11 is rather accurate with re- 
spect to the exact Onsager value T c ss 2.27. In addition, 



the transition exhibits a vertical jump to the disorder 
phase reminiscent of the rather abrupt vanishing of the 
corresponding Onsager second order transition. 

Accordingly, although the "dis/order" phase produced 
by the GUF - Metropolis combination is not physical, it 
is an intriguing result worth to be understood. 

In parallel, contrary to what could have been expected, 
combining Glauber dynamics to GUF restores the Ising 
single transition at T c as 3.09, which arises the question 
to determine why Glauber and Metropolis dynamics lead 
to different equilibrium sates when combined with GUF 
in the case of the 2-d Ising model? 

The various behaviors obtained respectively with GUF 
- Metropolis, GUF - Glauber, Onsager and Mean Field 
are shown together in Figure [51 which provides some co- 
herent picture about the effect of various approximations 
on departing from the exact treatment of the 2-d Ising 
model. At this stage, GUF needs more investigation to 
understand its nature and find the physical origins of 
those discrepancies with both the classical Mean Field 
and the exact result. The question has been evoked in a 
series of works 

One promising direction is certainly to extend the clus- 
ter size to which GUF has been applied. Instead of a 
cluster of 5 spins limited to the nearest neighbors, it 
should be interesting to consider both, the inclusion of 
next-nearest neighbors with a 9 spins cluster and also the 
extension to 3-d with a 7 spins cluster, noticing that the 
respective topologies are instrumental in the calculation 
of the various GUF coefficients. 
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